Mon. Not. R. Astron. Soc. 000, (1997) Printed 5 February 2008 (MN M^jX style file vl.4) 



Spectroscopy of the roAp star at Cir — II. 

The bisector and equivalent-width of the Ha line 



oo 

On 



I. K. Baldry 1 , M. Viskum 2 , T. R. Bedding 1 , H. Kjeldsen 2 ' 3 , S. Frandsen 2 

1 Chatterton Astronomy Department, School of Physics, University of Sydney, NSW 2006, Australia 
2 Institute of Physics and Astronomy, University of Aarhus, DK-8000 Aarhus C, Denmark 
3 Theoretical Astrophysics Center, University of Aarhus, DK-8000 Aarhus C, Denmark 



Accepted 1998 September. Received 1998 September; in original form 1998 May. 



Oh- 

CO 
CN 



OO 
CN 
On 
O 
OO 

6 



ABSTRACT 

We present bisector measurements of the Ha line of the rapidly oscillating Ap 
(roAp) star, a Cir, obtained from dual-site observations with medium-dispersion spec- 
trographs. The velocity amplitude and phase of the principal pulsation mode vary sig- 
nificantly, depending on the height in the Ha line, including a phase reversal between 
the core and the wings of the line. This supports the theory, suggested in Paper I, 
of a radial pulsational node in the atmosphere of the star. Blending with metal lines 
partially affects the Ha bisector results but probably not enough to explain the phase 
reversal. 

We have also detected changes in the equivalent-width of the line during the 
pulsation, and measured the oscillatory signal as a function of wavelength across the 
Ha region. 
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1 INTRODUCTION 

Rapidly oscillating Ap (roAp) stars are a sub-class of chem- 
ically peculiar magnetic (Ap or CP2) stars. They pulsate 
with photometric amplitudes below 8mmag and periods in 
the range 5-15 minutes, a Cir (HR 5463, HD 128898) is 
the brightest of the known roAp stars (V = 3.19). Previous 
observations of this star in photometry have shown that it 
has one dominant pulsation mode with a period of 6.825 min 
(/ = 2442 fiYLz; Kurtz et al. 1994, hereafter KSMT). Medupe 
& Kurtz (1998) observed that the amplitude decreased with 
increasing wavelength, from 2.71 ± 0.18 mmag in Johnson U 
(3670A) to 0.41 ± 0.13 mmag in I (7970A). They proposed 
that the rapid decline of photometric amplitude with wave- 
length could be explained by a decrease in the temperature 
amplitude of the pulsation with atmospheric height. 

Baldry et al. (1998a, hereafter Paper I) showed that 
the velocity amplitude and phase of the principal pulsation 
mode in a Cir varied significantly from line to line. However, 
it was difficult to interpret the data because of blending 
effects. In particular, there was evidence for a velocity node 
in the atmosphere of a Cir but it was uncertain whether the 
node was horizontal or radial. In this paper we look at the 
Ha line in more detail using the same set of observations, 
taken during two weeks in 1996 May. These observations 



comprise 6366 intermediate-resolution spectra taken using 
the 74 inch (1.88-m) Telescope at Mt. Stromlo, Australia and 
the Danish 1.54-m Telescope at La Silla, Chile (see Section 2 
from Paper I for further details). 

We first examined how the Ha profile changed during 
the principal pulsation cycle, in terms of the bisector at 
different heights in the line (preliminary results were pre- 
sented by Baldry et al. 1998b). In this way, the effect of 
the velocity on the profile could be analysed. To quantify 
the temperature effect, we measured the equivalent-width 
(EW) amplitude of the principal mode in filters of varying 
width and measured pixel-by-pixel intensity changes across 
the Ha region of the spectrum. Finally, we defined an ob- 
servable quantity (related to the EW of Ha) which had a 
high signal-to-noise ratio for the principal pulsation. Using 
this observable, some of the weaker modes in a Cir were 
detected. 



1.1 General Properties of a Cir 

Before we discuss the pulsation of a Cir, we will review the 
general properties in the light of a recent spectral analysis 
(Kupka et al. 1996) and the Hipparcos parallax measure- 
ment (ESA 1997). The distance of 16.4 ± 0.2 pc (parallax 
of 61.0 ± 0.6 mas), combined with a bolometric correction 
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of -0.12 ± 0.02 (M bol = 4.64, Schmidt-Kaler 1982), gives 
Mboi = 2.00 ± 0.04 (L = 11.4 ± O.4L ). From this lumi- 
nosity and the temperature of T e ft = 7900 ± 200 K (derived 
by Kupka et al.), we obtain R = 1.81 ± 0.11-Rq (angular 
diameter of 1.03 ± 0.07 mas). Combining this radius with 
log g = 4.2 ± 0.15 (Kupka et al.) gives a mass of M = 
1.9 ± 0.6AfQ. The rotation period of 4.48 days (derived by 
KSMT) and the radius means that v rot = 20.4 ± 1.2 kms -1 . 
Using vsvai = 13 ± lkms -1 (Kupka et al.), the inclination 
of the rotation axis to the line of sight is then i — 40° ± 5° . 
This is an improvement on the estimate given by KSMT due 
to more recent results. 



1.2 The Oblique Pulsator Model 

In some roAp stars, including a Cir, the amplitude of the os- 
cillations has been observed to be modulated by the rotation 
period of the star. In particular, the amplitude is at a max- 
imum when the observed magnetic field strength is also at 
a maximum. This led to the oblique pulsator model (Kurtz 
1982), in which an roAp star pulsates non-radially with its 
pulsation axis aligned with the magnetic axis. The ampli- 
tude modulation comes from the inclination of the pulsation 
axis with respect to the rotation axis. Therefore different as- 
pects of a non-radial mode are viewed as the star rotates. 
The oblique pulsator model predicts that, in the Fourier do- 
main, a frequency of mode £ is split into 21+1 frequencies. 
The frequency splitting is exactly equal to the rotation fre- 
quency, and the relative amplitudes are determined by both 
the inclination of the rotation axis to our line of sight (i) 
and the angle between the rotation axis and the pulsation 
axis (/3). Most of the review papers mentioned in the Paper I 
discuss the oblique pulsator model, a refinement of which is 
given by Shibahashi & Takata (1993). 

The principal pulsation mode in a Cir is believed to 
be a pure oblique dipole mode (£ = 1). KSMT observed a 
21 percent full-range variation of the amplitude during the 
rotation period. This variation implies that tan i tan j3 = 
0.21 ± 0.01 (see Section 4.1 from KSMT). Using our new 
estimate of i (see Section 



1.1), we get /3 = 14° ± 3°. This 
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means that the inclination angle of the pulsation axis to 
the line of sight (a) varies between 26° and 54° during the 
rotation cycle. 
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Figure 1. The Ha line in a Cir. The dotted lines divide the 22 
contiguous sections. 



'meaning 2nd/3rd order polynomial', low_rej=2, high_rej=4, 
niterate=20/25, grow=0. 

Initially, a sample of about 1400 points from each spec- 
trum was used to make a least-squares polynomial fit - 
3rd order for the Stromlo data and 2nd order for the La 
Silla data. A higher order fit was used for the Stromlo data 
because the continuum shape was less stable. Next, points 
below 2a and above 4cr from the fit were excluded from the 
next fit. The lower cutoff was used to exclude some absorp- 
tion lines and the upper cutoff only excluded cosmic ray 
events. The fitting and excluding routine was repeated 20 
to 25 times or in most cases until no more points were ex- 
cluded. In the final fit, about 1000 points were included, 
with none being within 50A of the core of the Ha line. Not 
all absorption lines were excluded from the final fit since 
there is very little real continuum at the resolution of our 
data (~1.5A). The real continuum level in each fitted spec- 
trum was around 1.005-1.010, depending on the wavelength. 
Therefore, we have scaled the spectra so that the continuum 
level is ~1.00 in the Ha region of the spectrum. The level 
varied slightly from spectrum to spectrum, with a standard 
deviation of about 0.001. 



2 GENERAL DATA REDUCTIONS 

Extraction of spectra and continuum fitting were done using 
IRAF procedures (see Section 3.1 from Paper I). We made a 
number of types of measurements on the reduced spectra, in- 
cluding bisector line-shift measurements which are described 
in Section 2^2 . The time series analysis of these is explained 



in Section 2.2 and is also applicable to other measurements, 
which are described later in the paper. Since the continuum 
fit is more critical than for Paper I, we first describe this 
procedure in more detail. 



2.1 Continuum Fitting 

The IRAF procedure continuum was applied to each 
spectrum in our data, using the following parameters: 
sample=@list, naverage=l, function=legendre, order=3/4 



2.2 Bisector Measurements of Ha 

The Ha line in each spectrum was divided into 22 contiguous 
horizontal sections (see Fig. [j] and Cols. 1-4 of Table |l[) and 
two extra sections for checking. For each section, the average 
wavelength of each side of the absorption line was measured, 
and a bisector line-shift (average position of the two sides) 
was calculated. We have used a telluric O2 band as a Doppler 
shift reference (see Section 3.3 from Paper I for details). 

An alternative method was also tried in which a least- 
squares fit was made to the position of each side of the line. 
A template spectrum was used to define the shape of the 
line and was shifted from side to side until a best fit was 
obtained. This method produced similar results and noise 
levels to the method of calculating the average wavelength, 
but the computational time was longer. 
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Table 1. Amplitudes and phases of the principal pulsation mode for the bisector velocity at different heights in the Ha line. The line 
was divided into 22 non-overlapping and two extra horizontal sections. For the results in this table, the Stromlo and the La Silla data 
sets were combined. 
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"Boundary of the section in relative intensity. 

6 Approximate mean width of the line at the height of the section. 
c Amplitude measured at 2442.03 /iRz. 

d rms- noise estimated from amplitude spectrum using the regions 1100—2300 /iHz and 2600-4400 /iHz. 

e Phase measured at 2442.03 /iHz, with respect to a reference-point (to) at JD 2450215.07527, with the convention that a phase of 0° 

represents maximum of the observed variable. 

f Error in the phase is taken to be arcsin (rms- noise/amplitude). 



2.3 Time Series Analysis 

For each section of the Ha line, the above analysis produced 
a time series of 6366 bisector line-shift measurements (4900 
from Stromlo, 1466 from La Silla). Each time series was 
high-pass filtered and then cleaned for bad data points by 
removing any points lying outside ±6.5 times the median 
deviation. Typically, about 100 data points were removed. 
Next, a weighted least-squares sine-wave fitting routine was 
applied (using heliocentric time) to produce amplitude spec- 
tra. 

For the analysis of the principal pulsation mode in 
a Cir, we have measured the amplitude and phase of each 
time series at 2442.03 ^tHz and estimated the rms-noise level 
by averaging over surrounding frequencies, 1100-2300 ^iHz 
and 2600-4400 ^iHz. The phases are measured at a temporal 
phase reference point (to) with the convention that a phase 
of 0° represents maximum of the observable, to was chosen 
to coincide with maximum light using data supplied by Don 
Kurtz (private communication); see Paper I for details. For 
the phase error, we have used simple complex arithmetic to 



find the maximum change in phase that an rms-noise vector 
could induce, i.e., arcsin (rms-noise/amplitude). 



3 BISECTOR VELOCITIES 
3.1 Results 

We have assumed that the bisector line-shift measurements 
represent velocities in the star. The velocity amplitudes and 
phases of the principal mode at different heights in the Ha 
line are shown in Fig. ^ and Table |l|. The results describe the 
oscillations in the bisector about the mean position at each 
height (see Fig. |^ for the bisector shape). From 0.4 to 0.8 
in the line, the amplitude decreases from 300 ms -1 to zero 
and then increases again, with a change in phase of —140°. 
Note that there is good agreement between the Stromlo and 
La Silla data sets, which were analysed separately for Fig. [| 
This gives us confidence that the observed variations are 
intrinsic to the star. 

In Paper I, the Ha velocity amplitude and phase 
were measured to be ~170ms _1 and ~330° using a cross- 
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Figure 2. Amplitudes and phases of the principal pulsation mode 
for the bisector velocity at different heights in the Ho line. Points 
with solid lines represent the Stromlo data and points with dot- 
ted lines represent the La Silla data. For each measurement, the 
vertical line is an error-bar while the horizontal line shows the 
extent of the section in the Ha line. 



correlation method. This is compatible with our bisector ve- 
locity results since the cross-correlation measurement would 
be dominated by the steepest part of the Ha profile, near 
the core. 

To further illustrate the behaviour of the Ha bisector, 
in Fig. ^ we show the bisector velocity as a function of time 
and height. The horizontal axis covers two cycles of the prin- 
cipal pulsation mode (P — 6.825 min). Only sections 0-19 
are shown (see Table |l]), since above a height of 0.8 the 
measur emen ts are significantly affected by line blends (see 
Section |3.2[ ). The figure clearly illustrates that the higher 
sections are pulsating nearly in anti-phase with the lower 
sections, and that the velocities of the middle sections are 
close to zero. Since the bisector velocity reflects the velocity 
at different heights in the atmosphere, these results support 
the hypothesis, suggested in Paper I, of a radial node in the 
atmosphere of a Cir. 



Figure 3. The Ha line in a Cir. The crosses show the mean 
position of the bisector. The wavelength is measured with re- 
spect to the centre of the line. The dotted line shows the shape 
of trapezium filter no. 2, which is used for the equivalent-width 
measurements. It is plotted using a different vertical scale. 
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Figure 4. Velocity field diagram. The arrows represent bisec- 
tor velocity vectors, relative to the mean velocity of the star, at 
different times throughout two pulsation cycles of the principal 
mode (P = 6.825 min). The velocity vectors are calculated from 
the amplitudes and phases of sections 0-19 given in Table |l| 



© 1997 RAS, MNRAS 000, §4l| 



Spectroscopy of a Cir - II. The Ha line 5 



3.2 Blending Considerations 

We attribute the velocity node (at ~0.65) and the phase 
jump between 0.4 and 0.8 to Ha line-formation effects, al- 
though there is some uncertainty as to how much blend- 
ing affects the results between heights 0.7 and 0.8. Above 
a height of 0.8, the results are significantly affected by line 
blending. In order to identify any metal lines in the wings 
of Ha, we used a synthetic spectrum calculated by Friedrich 
Kupka (private communication) and high-resolution spectra 
of a Cir taken in 1997 March (coude echelle, 74 inch Tele- 
scope, Mt. Stromlo). From this, we determined that most of 
the lines blended with the Ha profile from 6540A to 6585A 
are telluric, including all those below a height of 0.7 in the 
Ha line. 

At height 0.82, the anomalously large velocity ampli- 
tude is probably caused by the metal absorption line Sr I 
6550. 2A, with the possibility that Fe 1 6575. OA also con- 
tributes. Since metal lines can have amplitudes as large as 
1000 ms -1 (Paper I), this can explain the large bisector ve- 
locity seen at this height. At height 0.88, the amplitude 
is similar to that between 0.72 and 0.79 but the phase is 
anomalous. This measurement is affected by various metal 
and telluric lines, in particular a metal line blend including 
Mg II 6546.0A and Fe I 6546.2A. 

Below a height of 0.8, we expect any absorption fea- 
tures to have less impact on the bisector velocity due to the 
increasing steepness of the Ha profile. However, there is a 
metal line Fe I 6569. 2A in addition to a few telluric lines 
that affect the Ha profile between 0.7 and 0.8. This Fe I 
feature affects the bisector measurements between 0.71 and 
0.76. Therefore, it probably causes the jump in amplitude 
near height 0.72 shown in Fig. ^. It cannot explain why all 
the bisector velocities from heights 0.67-0.80 have phases 
around 190°-200°, but it does lower the significance of the 
phase jump between the core and this part of the wings. 
We speculate that without this Fe I feature, the Ha bisector 
velocity amplitude would increase steadily from 0.70 to 0.77. 

In summary, we believe features in Fig. below height 
0.7 reflect genuine effects in the Ha profile, while those above 
0.8 are dominated by metal and telluric lines. Between 0.7 
and 0.8, a blend probably exaggerates the phase reversal 
by causing an increase in the measured amplitude around 
height 0.73. 



3.3 Simulations 

Two sources of noise for the velocity measurements, espe- 
cially at high levels in the Ha line, are errors in the con- 
tinuum fit and changes in the total equivalent-width (EW) 
of the line (which are expected due to temperature changes 
during the pulsation cycle) . These sources will affect the bi- 
sector velocity if the line is asymmetrical at that height. To 
investigate, we first simulated fluctuations in the EW of the 
line and measured the resulting pseudo-velocity amplitude, 
and then simulated fluctuations in the continuum level. We 
concluded that, while there could be systematic changes in 
the velocity amplitude greater than 150 ms" 1 above a height 
of 0.65, the phase of the systematic change would be oppo- 
site between height 0.72 and 0.76. This means that this ef- 
fect cannot explain the different phase of heights 0.72-0.79 
compared to heights 0.40-0.56. Perhaps, this systematic ef- 



fect could explain the large difference in velocity amplitude 
between height 0.70 and 0.72. 

We also considered the effect of a systematic change in 
the slope of the continuum across the Ha line. This could 
arise because the continuum fit is not entirely independent of 
the metal lines and, since there are significantly more metal 
lines on the blueward side of Ha, the slope of the continuum 
may vary during the pulsation cycle. We have simulated the 
effect on the velocity amplitude of a continuum variation 
having a slope change of 0.001 / 100A. We found that such 
a large slope change could produce a systematic change in 
the velocity amplitude greater than 300 ms -1 above a height 
of 0.75 and thus mimic a velocity phase reversal. However, 
the simulation fails to reproduce a velocity node at 0.65 and 
other features of Fig. ^. Furthermore, from calculations of 
the intensity variations in various regions of the spectrum, 
we do not expect any systematic continuum slope changes 
to be larger than 0.0001 / lOOA. 

Our simulations show that continuum level, continuum 
slope and Ha EW changes cannot account for the observed 
features of Fig. [| Therefore, we conclude that the node 
and phase reversal are caused by the velocity field of the 
star. Note that random variations in the continuum level 
and slope may be larger than tested for, but this will only 
effect the noise level in the amplitude spectra. 

3.4 Discussion 

Hatzes (1996) has simulated line-bisector variations for non- 
radial pulsations in slowly rotating stars. His simulations 
show that a bisector velocity phase reversal could occur in 
modes with I — m > 3 (see Fig. 3 from Hatzes 1996). In fact, 
the Ha line-bisector variations in a Cir could closely be de- 
scribed by his simulations with £ = m = 3 or 4. However, he 
does not take into account any changes in pulsational am- 
plitude with depth and the simulations were done for much 
narrower metal lines. Given that there is strong evidence 
for the £ = 1 oblique pulsation model (KSMT), we suggest 
that a change in pulsational amplitude with depth is a more 
likely explanation for the bisector variations. 

It is interesting to note that similar behaviour to that 
seen in a Cir has been seen in the Sun. Deubner, Waldschik 
& Steffens (1996) have observed a phase discontinuity in 
the solar 3-min oscillations using spectroscopy. In particu- 
lar, they measured phase differences between the velocity of 
the core of the NaD2 line and various positions in the wing 
(called V— V phase spectra). They discovered a 180° phase 
jump in the V— V spectra near a frequency of 7000 /iHz. 
They also observed a phase discontinuity in the V— I (Line 
Intensity) spectra at a similar frequency. 

The 3-min oscillations are thought to be formed in 
an atmospheric (or chromospheric) cavity. Deubner et al. 
(1996) suggested a model to explain the phase discontinu- 
ities, involving running acoustic waves and atmospheric os- 
cillation modes, one of which must have a velocity node in 
the observed range of heights in the atmosphere. 

In the solar model, the eigenfunction of a p-mode with 
£ = 1 and n = 25 has a radial node separation in the outer 
part of the Sun of ~0.3 percent of the radius (Christensen- 
Dalsgaard 1998). In a Cir, if we assume the oblique dipole 
pulsation model (I = 1) for the principal mode, with a 
large frequency separation Afo = 50 ^tHz (KSMT), then 
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the overtone value of the principal mode is n = 48, using 
the asymptotic theory of asteroseismology (e.g. Brown & 
Gilliland 1994). The radial node separation in a Cir is then 
expected to be ~0.15 percent of the radius of the star. This 
is equivalent to a distance of about 1900 km, assuming the 
radius of a Cir is about twice the radius of the Sun (see 
Section [l.l| ). This does not take in to account the differ- 
ence in density as a function of radius between a Cir and 
the Sun. Using model atmospheres (Friedrich Kupka, pri- 
vate communication) to estimate the sound speed (v) just 
above the photosphere, we obtain a radial node separation 
of about 1500km (v/2f). This is in good agreement with 
the estimate obtained by scaling from a solar model. 

In the model of a Cir used by Medupe & Kurtz (1998), 
the geometric height between the formation of the I band 
and the B band is 250 km. We speculate that the extent 
of the line forming region is about 1000 km and that we 
are seeing one velocity node in the atmosphere. Recently 
Gautschy, Saio & Harzenmoser (1998) described pulsation 
models for roAp stars which suggest that radial nodes can 
be expected in the atmospheres of these stars. 
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Figure 5. The EW amplitude in filters of different width centred 
on the Ha line. The stars with solid lines represent the Stromlo 
data and diamonds with dotted lines represent the La Silla data. 



4 EQUIVALENT- WIDTH MEASUREMENTS 
OF Ha 

Kjeldsen et al. (1995) developed a new technique for detect- 
ing stellar oscillations through their effect on the equivalent- 
width (EW) of Balmer lines. We expect to find changes in 
the Ha EW in a Cir due to temperature changes in the star. 



4.1 Reductions 

We measured the EW changes of the Ha line directly by 
looking at intensity changes in regions of varying width 
across Ha. First of all, the spectra were linearly re-binned 
by a factor of 40 and shifted so that the centre of the Ha 
line was in the same wavelength position for each spectrum. 
This was to reduce the noise caused by instrumental shifts 
of the spectra. The value used to define the centre of the Ha 
line in each spectrum was taken from the cross-correlation 
measurements used for Paper I. Next, the mean intensities 
in 11 trapezium-shaped filters (Cols. 1-2 of Table were 
measured in each spectrum. Each filter was centred on Ha 
and the sloping part of the trapezium was two pixels wide 
(~ lA) at each end (see Fig. |^ for an example). 

For each filter, we obtained a time series of intensity 
measur emen ts which was analysed in the same way as in 
Section |2.3| to yield the amplitudes shown in Cols. 4-5 of 
Table From an intensity change (SI), we defined the frac- 
tional EW change as: 

SW _ -SI 
W ~ G-I 

where W is the equivalent-width, C is the continuum level 
(approx. 1.00 in our reduced spectra) and I is the mean in- 
tensity in the filter averaged over all the spectra (Col. 3 of 
Table ^). Using this formula, we converted intensity ampli- 
tude spectra to EW amplitude spectra. To test whether any 
bisector variations could affect these measurements, we sim- 
ulated shifts of ±600 ms -1 in the template spectrum. From 



this, we determined that the effect of such a signal on the 
filter intensities would be less than the rms-noise level. 

4.2 Results 

The results for the principal mode are shown in Cols. 6-10 of 
Table ^ and in Fig. |B| The phases of all the measurements lie 
between —10° and 15°. Since our phase reference point (to) 
coincides with maximum light (see Section 4 from Paper I), 
we conclude that the EW of the Ha line is pulsating in 
phase with the luminosity. This is expected, since maximum 
EW of the Ha line indicates maximum temperature in the 
stellar atmosphere. If the EW amplitude were the same in 
all filters then the change in intensity at each wavelength 
would be to be proportional to the depth of the absorption 
at that wavelength, i.e., 81 oc C — I . Such a profile variation 
is shown greatly exaggerated in Fig. [| Our results in Fig. ^ 
suggest that this is nearly the case, with an amplitude of 
1000 ppm, except that the core of the line is fluctuating in 
intensity by more than expected. 

We must remember that the measurements were made 
on continuum-fitted spectra, which means that the absolute 
fluxQ in the core might be constant while the wings and con- 
tinuum are changing in flux. We do have information about 
the pulsation in continuum flux. Medupe & Kurtz (1998) 
measured the photometric amplitude in Johnson R (central 
wavelength 6380A) to be 0.54±0.15 mmag (500 ± 140 ppm). 
A continuum flux amplitude of 500 ppm could only account 
for a 200 ppm Ha core relative intensity amplitude in contin- 
uum fitted spectra. This value is significantly smaller than 
the measured value of ~1000ppm (filter 1). Therefore, we 
can say with some certainty that the continuum flux from 
the star is varying in anti-phase with the flux in the core 
of the line and that, at some point in the wings of Ha, the 
flux is constant. To investigate this, we have calculated the 

* We use the term intensity when referring to continuum-fitted 
spectra and the term flux when considering the true flux from 
the star. 
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Table 2. Equivalent-width (EW) amplitudes and phases of the principal mode for filters of different width across Ha. T he m easured 
relative intensity amplitudes (SI) are converted to fractional EW amplitudes (SW/W), using the formula given in Section 1.1 For the 
results in this table, the Stromlo and the La Silla data sets were combined. 



filter 


filter 


mean 


intensity 


intensity 


EW 


EW 


S/N 


EW 


phase 


no. 


width 


intensity 6 


ampl. c 


noise 1 * 


ampl. c 


noise d 




phase 6 


error f 




(A) 




(ppm) 


(ppm) 


(ppm) 


(ppm) 




(°) 


(°) 


Q 


1.0 


0.378 


1085 


72 


1743 


116 


15.1 


354 


4 


1 


2.0 


0.402 


996 


42 


1667 


69 


24.0 


354 


2 


2 


5.9 


0.511 


692 


20 


1415 


40 


35.4 


357 


2 


3 


13.7 


0.616 


111 


20 


1148 


51 


22.6 


1 


3 


1 


17.6 


0.650 


390 


20 


1114 


56 


19.7 


3 


3 


5 


24.5 


0.694 


301 


20 


982 


65 


15.0 


1 


4 


6 


31.9 


0.730 


246 


20 


913 


76 


12.1 


7 


5 


7 


45.1 


0.779 


201 


21 


909 


95 


9.6 


9 


6 


8 


66.6 


0.829 


156 


21 


911 


122 


7.5 


13 


8 


9 


81.8 


0.853 


147 


21 


995 


141 


7.1 


12 


8 


10 


110.7 


0.882 


122 


20 


1039 


172 


6.0 


11 


10 



"Full-width half-maximum of filter centred on Ha. 

'Mean relative intensity in the filter averaged over all the spectra. The continuum level is ~1.00. 



c,d,e,f See TaWe |JJ 




plitude, we use the scaling law from Bedding et al. (1996) 
which can be written as; 



6540 6560 6580 

Wavelength (A) 



6600 



Figure 6. The Ha line with the dotted lines showing the profile 
with an increase and a decrease of 10 percent in EW. A variation 
of this type, but with much smaller amplitude (1000 ppm), was 
used to generate the dashed line in Fig. fei and would produce a 
constant amplitude in Fig. |i[ 



intensity amplitude of the principal oscillation m ode at each 
pixel in the spectrum, as we describe in Section 5.1. 



8W _ 1 SlnW (5L\ 



Model atmospheres can be used to estimate the change in 
EW of the Balmer lines as a function of temperature. From 
Kurucz's models of the Ha profile with log<? = 4.0 (see 
Table 8A from Kurucz 1979), SlnW/SlnT ~ 2.8 between 
7500 K and 8000 K. Using this value, we predict an EW am- 
plitude of SW/W = 1.0 x 10" 3 . 

The measured amplitude of 1000 ± 140 ppm (filter 9, see 
Table ^ is in excellent agreement with the prediction. How- 
ever, to some extent the agreement is fortuitous because in 
reality, neither (8L/L)hoi nor SlnW/SlnT are particularly 
well known. In the first case, since the photometric ampli- 
tude of a Cir varies with wavelength by more than expected 
for a blackbody (Medupe & Kurtz 1998), this will affect 
the bolometric luminosity amplitude. In the second case, 
the value depends significantly on the effective temperature 
and the accuracy of the model, giving a possible range in 
5 In W / S In T of 2-4. Additionally, there are uncertainties in 
the sensitivity of the EW amplitude to different modes (Bed- 
ding et al. 1996) . Finally, note that the velocity amplitude is 
predicted to be 120 ms -1 using the same scaling laws (Kjeld- 
sen & Bedding 1995), while the measured amplitudes vary 
between and 1000 ms" 1 (Paper I). 



4.3 The Equivalent-width Amplitude 

Using the photometric amplitude for a Cir of ~1.7mmag 
(Johnson B) measured by Don Kurtz (private communica- 
tion) around the time of our observations, we can calcu- 
late the expected EW amplitude of the principal pulsation 
mode. The photometric amplitude converts to (5L/L)boi = 
1.5 x 10~ 3 using the relation of Kjeldsen & Bedding (1995) 
with T G fT = 7900 K. In order to convert this to an EW am- 



5 PIXEL-BY-PIXEL INTENSITY 
MEASUREMENTS 

Another way of examining changes in the Ha profile is to 
look at relative intensity changes across the Ha region as a 
function of wavelength (pixel by pixel). 
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5.1 Reductions 

We coul d apply the time series analysis, as described in Sec- 
tion 2.3, to the intensity of each pixel in our spectrum. How- 



tude of 1000 ppm. This is the amplitude (cosine compon ent) 



ever, this is computationally intensive and, since we are only 
interested in the amplitude of the principal mode, is unnec- 
essary. Instead, the spectra were phase-binned at the prin- 
cipal frequency. The phase-binning was done so that Fast 
Fourier Transforms (FFTs) could be applied to the inten- 
sity variations, and therefore save calculation time for the 
measurement of the amplitude of the principal mode. 

A sample of 4718 spectra from the Stromlo data set were 
used to produce the phase-binned spectra. Spectra which 
prod uce d bad data points from the EW measurements (Sec- 
tion [4 . l| ) were not included. Each spectrum was linearly re- 
binned by a factor of 10 and shifted so that the centre of the 
Ha line was in the same wavelength position for each spec- 
trum (determined by cro ss-c orrelation) . This was similar to 
the reduction in Section 



4.1 



Additionally, the intensities in 
each spectrum were slightly adjusted, in order to remove 
low-frequency variations of the mean intensity across the 
Ha region (filter 9, 6522A-6604A). This was necessary m 
order to avoid introducing noise at the principal frequency 
due to the low-frequency intensity variations. Using this ad- 
justment, we obtained consistent results between the filter 
measurements on the individual spectra and on the phase- 
binned spectra. The spectra were phase-binned at the prin- 
cipal frequency using a weighted mean in each phase-bin (25 
phase-bins were used but this number was not critical). Fi- 
nally, a FFT was applied to the phase-binned series of inten- 
sities, to measure the amplitude and phase of the principal 
mode and its harmonics for each pixel. 



5.2 Results 

The intensity amplitudes of the principal mode and harmon- 
ics, as a function of wavelength, are shown in Fig. |^. Note 
that a signal is expected in the first harmonic with an ampli- 
tude of about 8 percent of the principal mode (KSMT) , with 
higher harmonics containing negligible signal. Therefore, the 
harmonics give a good indication of the noise level (the rms- 
noise is ~ 60 ppm) . From Fig. it is clear that the simple 
description of the profile variation given in Section 4.2 is 



inadequate and that the metal lines have an impact on the 
intensity variations. For instance, the peak blueward of Ha 
could be caused by the Sr I 6550. 2 A line, while the Fe I 
6569. 2A line may also have an effect (this is more evident 
in the next figure) . The signal from the core of Ha is strong 
while the signal from the wings is much weaker. 

To show the phase information, we have chosen to plot 
the cosine and sine components of the principal mode sep- 
arately (see Fig. |^). The phase of most pixels is near 0° 
or 180° and therefore most of the information is contained 
within the cosine component, which represents the relative 
intensity variations that are in phase or in anti-phase with 
the photometric pulsation. The amplitudes of this compo- 
nent are mostly negative which means that at maximum 
light, the relative intensity is at a minimum i.e. the EW of 
Ha is at a maximum. 

Two different simulations are also plotted: 

(i) The dashed line represents the relative intensity am- 
plitude assuming 81 = —A(C~I), where A is the EW ampli- 



expected from the profile variation described in Section 4.2 
(ii) The dash-and-dotted line represents SI — —AI/C 
where A is the continuum flux amplitude of 500 ppm 
(Medupe & Kurtz 1998). This line shows how the relative 
intensity amplitude (cosine component) would behave if the 
absolute flux in this region were constant (6520 A - 6605 A), 
and the continuum flux outside this region were pulsating at 
the measured value of 500 ppm. This is not a realistic simu- 
lation but we use it to estimate how the absolute flux in this 
region behaves. Where the solid line is above the dash-and- 
dotted line, the absolute flux is pulsating in phase with the 
continuum. We see that the absolute flux of the inner 10A of 
the Ha line is pulsating in anti-phase with the continuum. 



5.3 Discussion 

To evaluate these results, we have to consider systematic er- 
rors. Unlike for the bisector velocities, any systematic errors 
in the continuum level directly affect the EW and relative 
intensity measurements. The direction of any error would be 
the same at all wavelengths, equivalent to an offset in the 
cosine amplitude shown in Fig. ^. From looking at pseudo- 
continuum regions outside Ha, we estimate the intensity am- 
plitude of the continuum level to be 100 ppm in the contin- 
uum fitted spectra (cosine component of —100 ppm). This 
slightly changes the details of the Ha profile variation. 

Other systematic errors may be caused by residual 
Doppler shift signals that were not removed by the reduction 
process. It was not possible to remove such signals com- 
pletely because of the variation in velocity amplitude and 
phase between different lines and between different heights 
in the Ha line. This could account for some of the sine com- 
ponent of the amplitude. 

Ronan, Harvey & Duvall (1991) measured the oscilla- 
tory signal in the Sun as a function of optical wavelength. 
One of their wavelength regions included the Balmer lines 
H7 and H<5. More recently, Keller et al. (1998) looked at 
the H/3 line region in the Sun, using a similar technique. 
Both groups showed that the absolute flux oscillations in 
the Balmer lines were reduced to about 70 percent of the 
continuum signal. This differs sharply from a Cir, where 
the Ha core absolute flux is pulsating in anti-phase with 
the continuum. This may be due to a pulsational tempera- 
ture node in the atmosphere of a Cir, as was indicated by 
the results of Medupe & Kurtz (1998), and by the velocity 
node. The relationship between the temperature node and 
the velocity node is unclear with the present data. 



6 DETECTION OF OTHER FREQUENCIES 

In order to detect weaker modes in a Cir, we need an observ- 
able which has a high signal-to-noise ratio in the amplitude 
spectrum. Of the observables discussed so far, the one with 
the highest S/N for the principal mode is filter 2 from the 
EW measurements (S/N= 35, see Col. 8 of Table ||). How- 
ever, for the measurements of intensity in different filters 
across Ha, the noise level is approximately the same for fil- 
ters 2-10 (20 ppm, see Col. 5 of Table |§). This means that 
the noise is caused by variations in the continuum level and 
not by photon noise. Therefore, we can improve the S/N 
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Figure 7. The solid line represents the relative intensity amplitude of the principal mode, while the other lines show three higher 
frequencies (harmonics). 
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Figure 8. The solid line represents the relative intensity cosine amplitude of the principal mode, this component represents intensity 
changes in phase or in anti-phase with the photometric pulsation. The dotted line represents the sine amplitude. The dashed line is a 
theoretical amplitude with SI = — 10 _3 (C — and the dash-and-dotted line shows SI = —5 X 10~ 4 //C which represents approximately 
where the absolute flux amplitude would be zero (see text for details). 



for the principal mode by dividing the intensity in one fil- 
ter by another and thereby reducing errors caused by the 
continuum fit. The highest S/N (= 5f) was obtained by di- 
viding the intensity in filter 2 (FWHM ~ 6A) by filter 7 
(FWHM ~ 45 A). We call this observable R cw (ratio of Ha 
core to wing intensity). This is analogous to narrow / wide 
H/3 photometry. 

We produced a time series of these R cw measurements 



which was analysed in the same way as in Section 2.5 . The 
amplitude (5R CW ) of the principal mode was measured to be 
727 ppm, with an rms-noise level in the amplitude spectrum 
of 14 ppm (see Fig. []) . To search for other frequencies which 
were detected by KSMT (the same numbering is used), the 

© 1997 RAS, MNRAS 000, 



principal pulsation mode was then subtracted from the time 
series to produce a pre-whitened amplitude spectrum (see 
Figs |l 0| |l 1[) - In this spectrum, we have detected the modes 
fn, /s and If i, and the rotational splitting of the principal 
mode (/i ±r), with amplitudes greater than 2.3 x rms-noise. 
We have not measured the frequencies but have used the 
values given by KSMT, except we have increased the fre- 
quencies of /i to fs by 0.03 /iHz for consistency with the 
change in the principal frequency, as measured during our 
observations (see Paper I). The rotational splitting r is taken 
as 2.59 /iHz, as measured by KSMT. 

The amplitudes and phases of each mode are shown in 
Table ^[ with the last two columns showing the amplitude 
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Figure 9. Amplitude (SR C w) spectrum of the intensity ratio measurements. 




Figure 10. Amplitude spectrum after subtracting the principal frequency. The dotted lines represent the frequencies of modes detected 
by KSMT. 
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Figure 11. Same as Fig. [l^, at a higher frequency. The dotted line represents the frequency of the first harmonic of the principal mode. 
The inset shows the spectral window on the same frequency scale. 



© 1997 RAS, MNRAS 000, 



Spectroscopy of a Cir - II. The Ha line 11 

Table 3. The amplitudes (8R CW ) and phases of different modes measured using the intensity ratio of filter 2 and filter 7 (see Table 0). 
Except for the principal mode (/i), the measurements are made after pre- whitening with the principal frequency. 



mode 


freq. a 


ampl. 


S/N 6 


phase c 


phase 


relative amplitude 




\jinz) 


(ppm) 




1°') 




SR cw e 


photometric-^ 


h 


2265.46 


26 


1.8 


279 


3.3 


0.036 


0.055 


h 


2341.82 


20 


1.4 


104 


45 


0.028 


0.063 


h 


2366.97 


35 


2.4 


272 


24 


0.048 


0.057 


h - r 


2439.44 


37 


2.6 


92 


23 


0.050 


0.095 


h 


2442.03 


727 


50.8 


174 


1 


1.000 


1.000 


fi+r 


2444.62 


48 


3.4 


306 


17 


0.066 


0.115 


h 


2566.52 


34 


2.4 


16 


25 


0.047 


0.046 


2/i 


4884.06 


48 


3.4 


76 


17 


0.066 


0.077 



"Frequencies taken from KSMT. See text for details. 

'rms-noise estimated to be 14 ppm from the pre-whitened amplitude spectrum in the region 1100— 4400 /^Hz. 

c Phase measured with respect to a reference-point (to) at JD 2450215.07527, with the convention that a phase of 0° represents 

maximum of the observed variable. 

d Error in the phase is taken to be arcsin (rms-noise/amplitude). 
e Error of approximately 0.020 (rms-noise/amplitude of principal mode). 
^Stromgren v photometry (KSMT), with an error of approximately 0.007. 



relative to that of the principal mode, both for our spectral 
data and for KSMT's photometric data. A possible signif- 
icant difference is for the rotational side-lobes. Both side- 
lobes have a lower relative amplitude by more than 0.04 in 
our data. This rotational splitting is caused by a variation 
in the amplitude of the principal mode during the rotation 
cycle (P ro t = 4.479 days). Our results therefore suggest that 
the amplitude 8R CW of the principal mode varies by less dur- 
ing the rotational cycle than the photometric amplitude. 

To examine this, we divided the time series into 32 
shorter time periods of about 4 hours each. In Fig. [l^, we 
show the amplitude SR CW of the principal mode for each 
time period. The solid line shows the amplitude variation 
expected from the measurements of /i, fx — r and fi + r 
in the amplitude spectrum of the complete time series (Ta- 
ble ^) . The dotted line shows the expected amplitude varia- 
tion based on the full range variation of 21 percent and the 
ephemeris from the photometric results (KSMT). There is 
good agreement between the amplitude maximum of our fit 
and the ephemeris of KSMT, which supports the accuracy 
of the rotation period derived by them. Both fits shown in 
Fig. [li] appear to be consistent with the results of the anal- 
ysis using the shorter time periods. Therefore, we cannot 
say with certainty that the percentage amplitude variation 
is less in our data than the variation measured by KSMT. 
Such a difference, if it exists, may be due to limb-darkening. 

The average value of R cm itself in each time period var- 
ied slightly during the observation period and was slightly 
different between the Stromlo and La Silla data sets (see 
Fig. ^). The difference in the value between the two data 
sets was about 0.4 percent and was due to the inexact match- 
ing of the filters. The variation with time seen in Fig. [[3] is 
probably caused by the rotation of the star, which is slightly 
inhomogeneous in surface brightness. This variation does not 
significantly affect the measurement of 8R CW at pulsation 
frequencies. 

Variation of Balmer line profiles with rotation is com- 
mon in Ap stars. For example, Musielok & Madej (1988) 



investigated 22 Ap stars of which 17 showed variation of the 
H/3 index with a typical amplitude of 0.02 mag (19000 ppm). 
Our observed 5R CW amplitude of about 1500 ppm at the ro- 
tation frequency in a Cir corresponds to a much smaller 
Balmer line variation than is typical of Ap stars, even if we 
allow for some difference between the H/3 index and R cw . 
To test this, we calculated an intensity ratio using filters 
across Ha of similar widths to those used for the H/3 index 
(FWHMs ~ 40A and ~ 150A). The rotational amplitude of 
this intensity ratio was 1000-1500 ppm, which is more than 
a factor of ten smaller than the 0.02 mag found by Musielok 
& Madej in other Ap stars. 



7 CONCLUSIONS 

The existence of a radial standing wave node of the princi- 
pal pulsation mode has been put on a firm footing by the 
velocity phase reversal in the Ha line. However, there is only 
a 140° velocity phase shift between the Ha core and higher 
in the line. If we were seeing a pure standing wave in the at- 
mosphere then we would expect a 180° phase shift. A trav- 
elling wave compone nt, blending or the systematic effects 
described in Section 3.3 may cause the discrepancy. Addi- 



tionally, many of the metal lines studied in Paper I show 
phases that are neither in phase or in anti-phase with the 
majority of the lines. We proposed that some of the phases 
were anomalous and did not represent velocity phases due 
to blending effects. Therefore, in a time series of less blended 
(higher resolution) spectra, there should be a clearer distinc- 
tion between those lines that are formed above and below 
the velocity node and fewer lines with anomalous phases. 

Analysis of the Ha profile during the principal pulsation 
cycle, shows a large change in relative intensity of the core 
of the line. This means that the absolute flux of the core is 
pulsating in anti-phase with the continuum, indicating a pul- 
sational temperature node in the atmosphere, as suggested 
by Medupe & Kurtz (1998) on the basis of multi-colour pho- 
tometry. However, Medupe, Christensen-Dalsgaard & Kurtz 
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Figure 12. Amplitude (8R CW ) of the principal pulsation mode 
during separate time periods of approx. 4 hours. The asterisks 
represent the Stromlo data and the diamonds represent the La 
Silla data. The error bar for each amplitude measurement is the 
rms-noise level. See text for an explanation of the solid and dotted 
lines. 



homogeneities in the distribution of certain elements, and to 
test non-adiabatic non-radial pulsation equations. 

By defining an observable which detects the relative in- 
tensity changes in the core of the Ha line, R C w, we detected 
some of the weaker modes in a Cir. This is probably the 
first detection, in roAp stars, of modes with photometric 
(B) amplitudes of less than 0.2 mmag, using a spectroscopic 
technique. The ratio between spectroscopic amplitudes and 
photometric amplitudes can be used to identify the i value 
of different oscillations modes, as has been done recently by 
Viskum et al. (1998) for the 5 Scuti star FG Vir. We do not 
have enough S /N to identify the modes in a Cir, but we show 
the potential for using this mode identification technique. 

a Cir offers one of the best chances for theoreticians to 
model a star where magnetic fields are important both for 
its evolution and pulsation, because of the detailed infor- 
mation that can be gained from its oscillation modes using 
spectroscopy and photometry. Additionally, there is already 
a well determined luminosity using Hipparcos and, when an- 
gular diameter measurements are made, there will also be 
direct measurements of the effective temperature and radius 
of the star. 
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Figure 13. Average intensity ratio (R C w) during separate time 
periods of approx. 4 hours. The asterisks represent the Stromlo 
data. The diamonds represent the La Silla data, which have been 
increased by 0.4 percent to match the Stromlo data. The dotted 
line shows a sine- wave with a period equal to the rotation period 
of the star. 



(1998) show that non-adiabatic effects can explain the pho- 
tometric results without the need for a node. Perhaps, non- 
adiabatic effects have a significant effect on the Ha profile 
changes in a Cir. 

The next stage in the analysis of the mode dynam- 
ics in the atmosphere could be the use of high resolution 
(R~50000) spectroscopy to study the velocities, bisectors 
and EW changes of unblended metal lines, combined with 
model atmospheres to calculate the formation height of the 
lines. This may allow us to decide whether there is travel- 
ling wave component, to relate the velocity changes to the 
temperature changes, to measure any surface or vertical in- 
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